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After almost half a century since the work of Anderson [Phys. Rev. 109, 1492 (1958)], at present 
there is no well established theoretical framework for understanding the dynamics of interacting 
particles in the presence of disorder. Here, we address this problem for interacting bosons near T = 0, 
a situation that has been realized in trapped atomic experiments with an optical speckle disorder. 
We develop a theoretical model for understanding the hydrodynamic transport of finite-size Bose- 
Einstein condensates through disorder potentials. The goal has been to set up a simple model that 
will retain all the richness of the system, yet provide analytic expressions, allowing deeper insight 
into the physical mechanism. Comparison of our theoretical predictions with the experimental data 
on large-amplitude dipole oscillations of a condensate in an optical-speckle disorder shows striking 
agreement. We are able to quantify various dissipative regimes of slow and fast damping. Our 
calculations provide a clear evidence of reduction in disorder strength due to interactions. The 
analytic treatment presented here allows us to predict the power law governing the interaction 
dependance of damping. The corresponding exponents are found to depend sensitively on the 
dimensionality and are in excellent agreement with experimental observations. Thus, the adeptness 
of our model, to correctly capture the essential physics of dissipation in such transport experiments, 
is established. 

PACS numbers: 67.85.De, 03.75.Kk, 05.60.Gg 



I. INTRODUCTION 

Bose-Einstein condensates generated in trapped atom 
experiments constitute mcsoscopic quantum objects that 
are well localized in space and exhibit quantum proper- 
ties such as interference and phase coherence. They rep- 
resent a highly controllable and potentially rich template 
for addressing fundamental questions related to quantum 
measurement and configuring precision measurement de- 
vices [lj ■ Here we study the unique problem of transport- 
ing such finite-sized BEC's on rough surfaces, a situation 
inevitably encountered with BEC's generated on chips. 
There the roughness originates from the fluctuations of 
the surface fields 0). 

In a broader context, transport of superfluids through 
disordered potentials has always been a topic of theoreti- 
cal intrigue ever since the first experimental observations 
concerning super-flow of Hclium-4 through porous media 
Q. Huang and Meng first proposed a hard-sphere Bose 
gas in a random potential as a model describing the static 
properties of this scenario [H . Even though mostly quali- 
tative, this work was very insightful and provided definite 
clues to the connection between boson localization and 
condensate depletion. Subsequently, using a quantum 
hydrodynamic approach, Giorgini et al. Q looked at the 
disorder induced phonon damping in a Bose supcrfluid. 
However, not much else was said about dynamics. In fact, 
after almost two decades, questions related to the trans- 
port of superfluids through disorder have hardly been 
touched mostly due to the lack of proper experimental 
setups. However, very recent, highly precise experiments 
on trapped atoms 6 8] have begun to expose puzzling 



new aspects of such transport through disorder. A sharp 
reduction in dissipation rate beyond the Landau critical 
velocity Q is observed. Further, the experiments also in- 
dicate a power-law dependance of the dissipation rate on 
the interaction strength. Finally, when we view these dis- 
ordered BEC experiments as well localized macroscopic 
quantum objects sliding on rough surfaces, it touches 
upon another interesting research area with open ques- 
tions related to the origins of non- linear friction laws ■ 

Trapped atomic systems constitute a promising new 
platform whereby the many-body physics of various 
disorder-induced phenomena may be understood. Con- 
trary to condensed matter, these atomic systems are in- 
trinsically pure and allow for controlled insertion of dis- 
order by, e.g., modulating the lattice potential fllj ]. or 
creating an optical speckle pattern [f| Il2| . Experimen- 
talists are even able to control, to a good extent, the 
correlation properties of these disorder potentials. 

The significance of this type of extreme control needs 
no further discussion if we simply note that, for the first 
time, experiments are able to provide explicit access to 
the localized wavefunction, thus allowing direct evidence 
of Anderson localization, an exquisite phenomenon re- 
sponsible for the complete disappearance of metallic con- 
ductivity fTT| - [T3j . Such spectacular developments, with 
no surprise, give compelling reasons to theoretically in- 
vestigate the disordered Bose system in the light of the 
new techniques of probing made available by ultra-cold 
atoms. Especially, with the possibility of tuning the 
inter-particle interactions via techniques of Feshbach res- 
onance [h4T |. ultra-cold atom experiments provide a sys- 
tematic way to investigate the role of interactions in the 
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localization phenomenon; an aspect that has generated 
much debate in the field. Numerous transport experi- 
ments can be performed, requiring new insights due to 
the finite spatial extent of these trapped-atom conden- 
sates. In addition, excitations in a finite BEC are funda- 
mentally different in symmetry from those in the transla- 
tionally invariant counterpart. Unfortunately, at present 
there is no well established theoretical framework for 
describing such non-equilibrium hydrodynamic scenar- 
ios. Experimental observations p-Q continue to demand 
a clear understanding of the physical mechanism. Re- 
cent numerical simulations based on the Gross-Pitaevskii 
equation [l5j . confirm some of these observations and 
provide some hints towards the possible physical mech- 
anism, but a more detailed theoretical model is missing. 
As will be shown below, our model provides a micro- 
scopic underpinning of the experimental and numerical 
results. This allows us to determine exponents govern- 
ing the power law dependence of the dissipation on the 
scattering length, giving a deeper understanding of the 
interplay between disorder and interactions. Also, as will 
be evident from our discussion, such exponents critically 
depend on dimensionality. 

In this paper, we present a theory for determining the 
transport properties of a trapped BEC flowing through 
a disordered potential, similar to that used in the exper- 
iments of Rcfs. [(| 0, [Hf. In the past, similar interest 
in transport of localized wavepackets, for example in mi- 
crowaves [l6| , has been limited to time scales of the order 
of the diffusion time r. While our proposed framework 
may also be applied to this regime as well, we intend to 
focus on the specific experimental situation discussed in 
Refs. [H, 0. Here, the dissipativc transport is achieved 
by generating a center-of-mass mode of the trapped con- 
densate and allowing it to propagate through an optical 
speckle potential (see Fig.[T|). The damping of this mode 
is measured as a function of interaction strength and the 
center-of-mass velocity. The strength of the disorder is 
such that the period of oscillations T <r. 



II. COLLECTIVE MODES (EXCITATIONS) OF 
A TRAPPED BEC 



A tightly confined BEC can be made to oscillate as a 
whole, for instance, by abruptly shifting the trap by a 
small distance. If the displacement is small, these exci- 
tations can be easily described by linearizing the Gross- 
Pitaevskii (GP) equation 



above GP equation, interactions play an important role 
in defining these excitations. 
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FIG. 1: (Color online) Sketch of the experimental scenario in 
question. The zigzag line corresponds to the speckle potential 
that is superimposed onto the smooth harmonic potential. 
The condensate (shown in blue) is shifted to one edge of the 
trap and released at t = 0, resulting in dipole oscillations. 

Among other types of collective oscillations exhibited 
by the trapped condensates, of special significance is the 
dipole mode. It corresponds to the center of mass motion 
of the condensate which due to the harmonic confinement 
has the same frequency as the trap and is unaltered by 
interactions. In general, for 3D confinement, there exists 
a whole set of low energy modes that can be obtained 
analytically by solving the coupled hydrodynamic equa- 
tions (derived from the GP equation above) for the veloc- 
ity, v(r,i), and the density, n(r, t): dn/dt + div(vn) = 
and mdv/dt + V {mv 2 /2 + V cxt + An} = 0. These equa- 
tions allow for solutions, where the condensate is oscil- 
lating as a whole [l7[, with frequencies given by u> e — 
<-<Jho\/2n2 + 2n r £ + 3n r + I, assuming T4xt(r) is harmonic 
with frequency uj^o- The dipole mode is the one with 
quantum numbers (n r = 0, £= 1). In the absence of any 
other potential (disordered or not), the harmonic trap 
is special, leaving the center-of-mass motion completely 
decoupled from the relative degrees of freedom even for 
large amplitudes. However this feature is immediately 
lost as soon as any other external potential, e.g. disor- 
der, is turned on, resulting in coupling of different modes, 
eventually leading to the complete decay of a well defined 
oscillating mode. 



III. KEY CONCEPTS 

Our model consists of two main ingredients, disorder 
and inhomogeneity, which we describe below. 
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where V cx t(j) is the trapping potential, and the inter- 
action parameter A, in terms of the s-wave atom-atom 
scattering length, a, is A ~ AnT^a/m. The set of Bogoli- 
ubov equations that follow provide the eigenenergies, A^, 
of the elementary excitations. As is apparent from the 



A. Disorder 

For treating disorder within an analytic ap- 
proach, among several well established techniques 
in the literature [l8l . [l9| . we find it convenient 
to use the replica method. The disorder poten- 
tial E/d(x) is described in terms of a Gaussian 
distribution function with a strength 7 such that 
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P[U d ] = cxp[-l/(2 7 2 )/dxdx'[/ d (x)A- 1 (x- x')C/ d (x')] 
and K describes the spatial correlations [l8j|. Such 
procedure automatically allows for a freedom in choos- 
ing the actual disorder potential, the details of which 
are quickly erased by multiple scattering. To further 
simplify our investigation, we consider a white noise 
correlation function A(x) ~ S(x). Thus, (J7d( x )}dis = 
and (J7 d (x)[/d(x'))di s = 7 2 <5(x-x'), where (■ • -) dis means 
averaging over disorder. However, the trick of using the 
replica method lies in following a particular procedure 
for taking the sequence of averages (quantum followed 
by disorder), essential for calculating observables. While 
the replica trick has been typically used in relation to 
Anderson localization of non-interacting electrons, we 
adopt it in the context of ultra-cold interacting atoms. 
A brief outline of this technique is provided in the 
Appendix. Aside from technical aspects, the replica 
trick allows for a systematic diagrammatic perturbation 
expansion within which all transport properties can be 
calculated, at least in principle. 
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FIG. 2: (Color online) Sketch showing the different relevant 
distances defined in the text. Dark (light) blue color repre- 
sent spatial regions of the condensate where the local sound 
velocity is greater (lesser) than the center-of-mass velocity. 
Graphic arrows indicate excitation moving along the z direc- 
tion. 



IV. THE MODEL 



B. Inhomogeneity 

Next we use local density approximation (LDA) [l?} 
to account for the inhomogencous density distribution of 
the trapped condensate. This is achieved by defining a 
position dependent chemical potential, 

fj, -> p{r) = fi- Vtrap(r), (1) 

clearly justified in the Thomas- Fermi limit [13], which is 
typically always the case for trapped-atom BEC experi- 
ments. Thus, the density of the stationary condensate in 
the absence of disorder is given by 

P (r) = M (r)/A (2) 

Further, we assume that as the condensate undergoes 
dipole oscillations, on average the shape of the conden- 
sate is retained, meaning p(r + r cm (t)) = p(r). This is 
an extremely good assumption, it stems from the rigid- 
ity associated with the Bosc condensate. This greatly 
simplifies the problem by allowing us to define all the 
relevant quantities locally, for example we define the lo- 
cal Bogoliubov excitation energies [2fJ 

A k [r] = yfel AV(r) 2 (3) 

with efe = h 2 k 2 /2m + Ap(r). Also, from the linear part 
of the above dispersion relation, we can define a local 
sound velocity 



c[r] = y/Xp(r)/m. 



The total Hamiltonian for the problem can be written 
as the sum 

H = -fffree + Hj n i +Vtrap + U<i 
= H e fi + Vt ra p + ^Ud^ 

LDA replica, 

where Hf rco is for the free non-interacting atoms and H ni t 
is due to the atom-atom interaction. As indicated above, 
we first write an effective Hamiltonian for the bosons 
in the absence of any potentials by treating H- lnt at the 
mean field level. Diagonalizing H c ff, the lowest eigen- 
value corresponds to the condensate mode at tempera- 
tures T < Tbec • The excitations above the ground state 
arc nothing but the Bogoliubov modes described above 
and are orthogonal to the condensate mode. Any other 
external potential will in general couple these different 
eigenmodes. However, as mentioned earlier, the trapping 
potential is typically smooth and LDA holds, drastically 
simplifying the calculation. The only piece that remains 
to be included is the potential, Ud- At any spatial lo- 
cation, the disorder tends to deplete the condensate by 
populating the Bogoliubov modes. However, this does 
not happen if the condensate is at rest, with zero center- 
of-mass velocity. Due to the superfluid property of the 
BEC, Landau criterion holds and there is a non-zero crit- 
ical velocity below which no excitations can be generated 
and the noncquilibrium state will persist forever. This is 
precisely where the finite size of the condensate plays an 
important role. At any given center-of-mass velocity v, 
on the contrary, there is always a finite region of the con- 
densate, shown by the light blue color in Fig. [3J where 
v > c[r] thus satisfying Landau criterion for excitations 
and providing a channel for the energy to be transferred 
from the center-of-mass motion to the Bogoliubov modes. 
In essence, this is the physical mechanism that provides 
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the frictional force responsible for slowing down the con- 
densate center-of-mass motion. 

Remark: Before we proceed to calculations, we need to 
take into account the fact that, previously, excitations 
were calculated in the rest frame of the condensate. Thus, 
if we wish to continue calculating everything else in the 
same frame, it implies having to use a time dependent 
disorder potential. 

This is certainly a non-trivial complication and cannot be 
simply neglected. Fortunately, LDA turns out to be per- 
fectly suited for such a situation. It allows us to work in 
the lab frame instead, hence static disorder. Now, from 
Eq. ([3]), the locally defined Bogoliubov modes are eigen- 
states of momentum and hence remain unaltered under 
a Galilean frame transformation, only their energies are 
shifted according to: 

A k ^A|f b = A k -Mc-v. 

Thus, we have gathered all the ingredients for conduct- 
ing the transport calculations via the diagrammatics of 
the replica trick. The above model is quite general and 
allows for calculating obscrvablcs in a broad range of sce- 
narios by including appropriate orders in the perturba- 
tion expansion. For example, if the time scale of interest 
is larger than the Heisenberg time (time scale defined by 
the trapping potential) , one needs to include appropriate 
particle-hole diagrams to describe the diffusive dynamics. 
On the other hand, if the dynamics occurs on the trap 
scale, then depending on the disorder strength, appropri- 
ate single particle diagrams may be sufficient. The only 
point one needs to remember is that the bare propaga- 
tors of the replica theory are those for the quasiparticlcs 
representing atoms dressed by interactions with the cor- 
responding Green's function given by [2(| [2l[ : 
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Sf(k,«j„)[r] - lab , lab 



1 



A{f[r]-iw n A{f[r]+W 
where ui n is the bosonic Matsubara frequency. 

V. EXPERIMENT 

As elaborated in the introduction, our model is moti- 
vated by recent experiments performed on cold trapped 
atomic condensates. Hence, we illustrate the applicabil- 
ity of the model by considering the actual experimental 
situation of Refs. [H, 0|: 7 Li atoms trapped in an ax- 
ial harmonic trap with u z = 2n x 5.5 Hz, aspect ratio 
a = lo±/oj z ~ 46. The s-wave scattering length is tuned 
via a Feshbach resonance to about a = 25 ae, resulting 
in an axial condensate with the size represented by the 
Thomas-Fermi (TF) radius, which in units of the axial os- 
cillator length £ z , is given by R sa 13.6 £ z . These param- 
eters conform with the experimental data for damping of 
the BEC dipole mode used in the plot of Fig. [3] The dis- 
order potential is produced using an optical speckle with 



gaussian correlation: 

([/ d (x)[/ d (x'))dis = V^exp 



2(x-x') 



i\2 



(4) 



Further, while no data on temperature is available, from 
the observations presented, it is safe to assume that 
the experiment is conducted at temperatures well below 
Tbec and thus aT = theoretical prescription should 
suffice. 

The TF density of the stationary condensate in cylin- 
drical coordinates follows from Eqs. ([1]) and @ 



P{r±,z) 



muJ z (R — a r\ 
2A 



with the normalization fixed by the TF axial size R via 



N = 4tt 



R/a rR[r ± ] 



p(t±, z)dz rj_dr±, 



(5) 



with the different distances indicated in Fig. [2j For 
convenience, we define the dimensionless ratio £[r±] = 
v/c[r±], where c[r±] is the speed of sound at the coor- 
dinate {r±,z = 0}. We also use the notation, c = c[0] 
and £ = £[0]. We point out that even though the aspect 
ratio is large, the chemical potential yu ~ 0(hu±), hence 
a full 3D treatment is required. Finally, we observe that 
the oscillation data shown in Fig. [3] indicates dynamics 
over time scales much shorter than the diffusion time 
dictated by the disorder. This greatly reduces the ma- 
nipulations and it suffices to calculate only the single par- 
ticle diagrams. The imaginary part of the selfenergy, ob- 
tained from the dressed propagator of the non-interacting 
quasiparticles, immediately provides the transition rate 
W(k — > k',r) at the spatial point r (for details refer to 
the Appendix). However, we alert the reader that one 
typically never measures the single particle propagator 
in an experiment and we still need to connect to the 
damping of the classical dipole mode. We achieve this 
by recalling that the energy per condensate particle is 
nothing but the chemical potential /i(r). Thus, the total 
energy transfer rate out of the k = condensate mode is 
precisely 



(6) 



r = / dr / dk fj,(r)W(0 -> k,r), 



where k is some infrared cutoff imposed by the trapping 
potential and Q is the condensate volume. The above in- 
tegral can be approximated and has a nice analytic form 
if excitations are allowed only along the axial (z) direc- 
tion. This constraint can be understood to arise from 
the high aspect ratio a of the trapping potential. Even 
at large fc, where this approximation may fail, we expect 
the density of states to be diminished enough to accrue 
any appreciable error. With this we arrive at the final 
result in terms of dimensionless quantities represented 
by (~)(physical quantities scaled in units of the axial 
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harmonic trap: lengths in units of I 
energies in huj z ) 



fR 4 



ttcl Jq 



\/Kpim>Jz and 
t 2 2[Z/Vi,Un/Vi}dt, (7) 



where the upper (lower) integration limit corresponds to 
the center (radial edge) of the condensate, £ m j n is related 
to k, £ m in = hn/2mc and the function X is defined by 
Eqs. dlT|) and flH2} of Appendix. 



VI. RESULTS AND DISCUSSION 

To begin with, we first look at the plot in Fig. 03a) 
that shows the velocity dependence of the function 1 for 
points along the long axis of the condensate. The sharp 
discontinuity near the speed of sound c is a clear indi- 
cation of the origin of distinct regimes of damping as 
the BEC flows through the disordered medium. Essen- 
tially, for velocities significantly smaller than the speed 
of sound c, excitations are only allowed in a small shell 
near the surface of the condensate. This region tends to 
grow inwards, eventually allowing excitations everywhere 
in the condensate, as the center-of-mass velocity reaches 
c. However, increasing the velocity beyond this critical 
value, the excited volume remains fixed while the Bogoli- 
ubov density of states keeps falling, resulting in net decay 
of the damping rate with velocity. 

The above expression for T represents the first of the 
two main results of this article. It conveniently leads us 
to a dynamical equation for the condensate motion, if we 
picture the BEC oscillating in the random potential as a 
rigid body whose hydrodynamic oscillations are damped 
as the kinetic energy is transferred to the internal ex- 
citations. Such damped oscillatory motion of the BEC 
center-of-mass coordinate can easily be described by an 
equation for the peak velocity of the condensate motion 
in terms of £ pca k 



peak 



at 



NR 2 ^ 



peak 



- -T[£pcak] ~~ /^Cpcak; 



(8) 



where the energy dissipation rate is averaged over one os- 
cillation period, represented by f = (1/T) J Q T Tdt. The 
above equation for £ pea k represents the second main re- 
sult of this article. These results immediately provide im- 
portant asymptotic exponents governing the dependence 
of the velocity-damping rate j3 on the scattering length. 
For this, all that is required is to remember that the to- 
tal number of atoms TV is fixed, automatically implying 
R ~ a 1 / 5 from the normalization condition Eq. (JSJ; 
thus, c ~ a 1 / 5 follows as well. Now, from the functional 
form given by Eq. (|B2|) . 



-6/5 



V£« 1 
V£> 1 
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(9) 



The latter limit is easily confirmed from the plot, Fig. 10, 
in Ref . Q . While, data for the former limit is not present 




FIG. 3: (Color online) (a) Damping function along the long 
axis of the condensate plotted as a function of f . (b) Thick 
solid curve: theoretical prediction for the peak center-of-mass 
displacement of the condensate as a function of time for 7 = 
0.5, filled circles connected with thin line are experimental 
data points for the center-of-mass position as the condensate 
oscillates in the trap Q. Square marker indicates the point 
where £ peak = 1. 



in the same figure, available data points for £ < 1 show a 
clear trend towards the predicted limiting behavior. In- 
stead, we indirectly verify this limit by extracting the 
damping rate as a function of £ from the data shown in 
Fig. (3](b) which can then be written as a function of scat- 
tering length, a, by fixing v in Eq. ([?])■ We again find very 
good agreement, confirming the exponent corresponding 
to the former limit as well. 

The only parameter that remains to be estimated is the 
5-correlated disorder strength 7. We remind the reader 
that, in the actual experiment, the disorder potential, 
Eq. (jlj, has a finite correlation length a. Such spatial 
correlations arc indeed outside the scope of the current 
LDA based model. Therefore, 7 in our model maybe con- 
sidered as a free parameter representing the renormalizcd 
value of the disorder strength. Our result, the solution 
of Eq. © plotted in Fig. 03b), is in excellent quanti- 
tative agreement with the experimental data, indicating 
the validity of a delta-correlated ansatz for the speckle 
potential. The proposed model is also able to capture all 
the damping time scales via only coupling to the Bogoli- 
ubov excitations. Thus, we identify the most dominant 
source of dissipation, sufficient to describe transport in 
such ultra-cold trapped- atom experiments. 

Finally, we shall comment on the strength of the dis- 
order. For this, let us use the value of 7 estimated in 
Fig. 0£b) to find the magnitude of disorder-induced con- 
densate depletion given via the method of Ref. [3| 



Pdep(O) 



2 2 
m 7 



8tt 3 / 2 ^ 4 V a 



P(0) 



We find pdcp(0)//o(0) = .002, implying that the conden- 
sate density remains almost unchanged after inclusion of 
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speckle potential and thus justifying the assumption of 
weak disorder, permitting perturbative treatment up to 
second order in disorder strength. Furthermore, it also 
points to an important aspect about the disorder poten- 
tial used in the experiment Q. There, even though the 
disorder strength seems to be large, given by Vd ~ 0.6 fj,, 
our theoretical model points to a critical interplay be- 
tween disorder correlations and interactions, resulting in 
a dressed disorder with a small rcnormalizcd strength 7. 



VII. CONCLUSION 

To summarize, we have given a theory for describing 
the transport properties of finite-size atomic condensates 
flowing through disorder potentials. We emphasize the 
intricate role played by the inhomogeneous character of 
the condensate, resulting in a completely different form of 
the damping function, X, above and below the critical ve- 
locity, c. The validity of our model is established via com- 
parison with the experimental observations of damped 
hydrodynamic dipole oscillations of the condensate Q- 
The predicted power-law dependance of the damping on 
the interaction parameter is in excellent agreement with 
experimental observations. Our results, while based on a 
white noise disorder model, show extremely good agree- 
ment with experimental data, thereby, indicating the in- 
significant role of finite correlations in determining the 
dissipation time scales. While a full theory was not 
necessary for the aspects considered here, an improved 
quantitative picture can be realized by including finite 
disorder-correlations, for example, via a coherent poten- 
tial approximation [22[ . Such procedure leads to an effec- 
tive, renormalized disorder vertex, whose value can then 
be compared with the 7 obtained here. However, these 
calculations are beyond the current model and would be 
considered in a future work. Furthermore, although the 
effects of Anderson localization turn out not to be rel- 
evant for these experiments [f| 0], we predict they could 
be if the center-of-mass motion is very slow and there- 
fore the physics is dominated by diffusion. These effects 
could be incorporated into our calculations by including 
appropriate diagrams corresponding to particle-hole pro- 
cesses. Finally, the analysis presented here can be easily 
translated in the language of friction indicating the man- 
ifestation of non-linear friction on mesoscopic quantum 
objects. 
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Appendix A: Replica Method 

The difficulty of disorder averaging is eliminated using 
the standard replica trick of considering R replicas of 
the same field with the understanding that the R 
limit is taken at the end. However, this procedure comes 
with a cost, it leads to an effective attractive interaction 
between atoms from different replicas. To illustrate this, 
we simply write down the full average of some operator 
6, as 



(0) 



lim 



1 d(Z R [r,]) dh 



R dr. 



(Al) 



where rj is some source field, and Z R [rj\ is the parti- 
tion function of the replicated boson field {(ip a<! , ip a ), a = 
1, .., R}. Thus, essentially it involves an averaging of the 
partition function over the distribution, P[{7d], that can 
be carried out quite easily using Gaussian integration re- 
sulting in an effective action 



R 

a,b=l 



with 



replica 



2_ 

' 2 



J2 / dx ^(^(^(^(4 



where S'ciean is the action in the absence of disorder, and 
index m and n are introduced to indicate Matsubara 
frequencies u> rn and uj n , respectively. Thus, the addi- 
tional replica-induced action represents interactions be- 
tween atoms of different replica index with the disorder 
vertex represented by the Feynman diagram of Fig. HJa) ■ 
Also, it does not involve energy exchange between repli- 
cas since the disorder potential is intrinsically time inde- 
pendent. 



Appendix B: Damping Function 

The relevant diagrams are the single particle ones as 
shown in Fig. EJb)-(c). A technical point of interest is 
that the contribution of diagrams containing loops of the 
type shown in Fig.|4jb) tend to zero due to the limit R — > 
in Eq. (|A1I) . In essence, the replica trick eliminates 
all diagrams that arc disconnected before the disorder 
averaging, thus vastly reducing the computation. 

For the particular experimental situation considered 
here, it suffices to include just the leading order diagrams. 
Further, the T = nature of the experiment makes the 
use of replica trick redundant and the result coincides 
with the Fermi's golden rule. However, we remind the 
reader that the model itself is more general and maybe 
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applied to other situations as mentioned in the text. All 
throughout, our goal has been to provide a theoretical 
model to describe transport in trapped, disordered, ultra- 
cold atomic experiments in general. Thus, the formal 
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FIG. 4: Diagram contributions to the self energy: (a) corre- 
sponds to the disorder vertex, m and n axe replica index, (b) 
and (c) are second order, (c) and (d) are fourth order con- 
tributions to the selfenergy. The external lines are not to be 
included for selfenergy calculations. The solid lines represent 
the quasiparticle propagator. 

expression for the transition rate is given simply by the 
imaginary part of selfenergy (diagram of Fig. H{c)). 

27TT 2 

W(0 -> k, r) = _J_p(r)(« k [r] - z; k [r]) 2 5(A k [r] - Hk ■ v), 



where Uk and v k are the amplitudes that define Bo- 
goliubov quasiparticlcs and are related to the energy 



A k = y/ei - X 2 p 2 by 



(e k /A k + l)/2 and v 



{e k /A k - l)/2 with e k = h 2 k 2 /2m + Xp. Therefore, the 
full expression for the loss rate T can be written simply 
by substituting in Eq. ([6]) of text. Performing the z and 
the k integration we arrive at 

2_-^2 2 3 rR/a 



where i?[r_i_], indicated in Fig. 2, is written as = 
y/ R 2 — ol 2 t 2 ^ and the function I is defined in terms of 
the boundary function T 



z\p, q) 

with 



F{l,p) - JVl -P 2 + Q 2 ,P), Yp < s/T+ q 



Vp> y/l + g 2 

.2 



(Bl) 



y(-5 + 2y 2 - 3a; 2 ) yj -\ + x 2 + y 2 + 



(3 + 2a; 2 + 3a; 4 ) log foy + 2^-1 + y 2 + x 2 ) 1(B2) 



Now, with straightforward manipulation, one can imme- 
diately show that Eq. ([7} of text follows by simply ex- 
pressing all quantities in dimensionless units. 
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